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15 CM THRUSTER PERFORMANCE TESTING 

Apparatus 

A SERT II thruster 1 was modified to include independent main and 
hollow cathode flow control and measurement systems, electromagnets, a 
solid disc baffle, and a cathode with a 0.05 cm dia. orifice. A larger 
cathode orifice than the one used on SERT II (0.015 cm dia.) was required 
to facilitate sustained operations at the higher beam current conditions 
expected in these tests. Power supplies and Langmuir probes used in the 
testing are described in Reference 2. Tests were conducted on the 
thruster described above with both standard flat grids and dished grids 
over a range of total flow rate conditions to facilitate performance 
comparisons. The two grid sets used in the tests are characterized by 
the dimensions shown in Figure 1. The dished grids, which are both 
0.046 cm thick, are seen to have 0.20 cm dia. aperatures in both the 
accel and screen grids on 0.25 cm centers, to be separated by 0.079 cm 
and to be dished out 0.5 cm at their centers. The flat grid spacing on 
the other hand is 0.23 cm and the holes in the 0.15 cm thick accel grid 
are slightly smaller than those in the .074 cm thick screen grid. The 
open area fractions for the screen grids are comparable. In spite of 
the greater voltages applied to the flat grids, their greater separation 
and hence lesser perveance effects an overriding reduction in the Childs' 
law beam current density to 32% of that calculated for the dished grids. 
The effective ion acceleration distance used in this calculation of 
Childs' law current densities is that suggested by Rawlin 3 for 30 cm 
dished grids. Since the beam current at the nominal operating condition 
decreased from 650 ma for the dished grids to 240 ma for the flat ones 
(a ratio of 38%) the Rawlin expression appears to be valid for this 
smaller thruster too. 

Performance 

The discharge losses (including keeper discharge losses) for the two 
grid systems are presented as a function of propellant utilization 
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efficiency in Figure 2 over a total mercury flow rate (m-j.) range from 
220 ma to 730 ma equivalent. The cathode flow rates (rh c ) , also indicated 
on this figure* were selected to produce an arc voltage in the range 
36-37 v (designated by the solid symbols) near the knee of each perform- 
ance curve. The data of Figure 2 obtained with the high perveance, 
dished grids show that the thruster exhibits a decrease in maximum utili- 
zation as the total flow rate is throttled from 730 ma through 440 ma to 
220 ma as Kaufman's model for maximum utilization 4 predicts. With the 
flat grids however the maximum utilization is greatest at the 440 ma 
total flow condition. A decrease in maximum utilization was observed 
with the flat grids as the flow rate was increased above about 440 ma 
several times, and it became so severe in some instances that beam current 
actually decreased as flow rate was increased at essentially constant arc 
voltage and current conditions. Comparison of the data at 440 ma total 
flow for the two grid systems indicates flat grids yield superior perform- 
ance at this lesser flow rate. This leads one to the conclusion that 
optimum thruster performance over a throttling range can be achieved by 
reducing the grid system transparency as flow rate is reduced. The trans- 
parency of a grid system is primarily a function of the screen grid open 
area fraction and the shape and location of the thruster plasma-beam 
plasma interface 5 , which are determined in turn by the effective ion 
acceleration distance. The screen grid open area fraction is not signifi- 
cantly different between the two grid sets, but the screen grid thickness 
and grid separation distance have been changed and they alter the effec- 
tive ion acceleration distance. This suggests optimum thruster perform- 
ance might be achieved by varying grid separation as the thruster is 
throttled, although tests have not been performed where grid separation 
alone is varied to achieve this optimum condition. The range of throt- 
tling suggested by the two sets of curves in Figure 2 is from a 175 ma 
beam current (80% utilization and 300 ev/ion) to 630 ma beam current 
(85% utilization and 250 ev/ion) - a 3.5 to 1 throttling range. 

The ion beam profiles measured with a Faraday probe swept through 
the beam at a location 5 cm downstream of the accel grid are shown in 
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Figures 3, 4 and 5. Integration of the current density profiles shown 
over the beam cross sectional area yields beam currents which agrees 
with the measured beam currents to within 5%. The parameter shown with 
each curve is the ratio of average to maximum beam current density appro- 
priate to that curve shown* Since an e near unity is desirable. Figure 3 
suggests the flat grids yield the better ion beam profile at the high 
flow rate. The parameter e is less for the dished grids because of the 
significant current densities at radii beyond 8 cm. This is a result of 
poor focus of some of the dished grid beamlets, due to grid hole misalign- 
ments - a condition that has since been corrected by redesign of the grids 
to compensate for the inherent misalignment which occurs when grids are 
dished. 

Comparison of Figures 3, 4 and 5 suggests the ion beam profile be- 
comes more peaked as the thruster is throttled and that this effect is 
most pronounced in the flat grid system. The tendency toward a more 
peaked profile as the thruster is throttled is considered to be due to 
the increased fraction of the propellant flow that is introduced through 
the cathode as it is throttled. This tends to cause an increase in 
neutral density near the center of the discharge chamber, hence greater 
ionization on the centerline and the more peaked profile. 

Plasma Property Contours 

Figures 6,7 and 8 present the electron temperature, plasma potential 
and primary and Maxwellian electron densities within the thruster over 
the range of flow rates investigated and at the 36-37 v operating condi- 
tions identified by the solid symbols on Figure 2. Each contour shows 
the variation of the indicated property on the surface of a plane passing 
from the thruster centerline to the anode over the length of the thruster. 
The cathode pole piece is identified by the line parallel to the thruster 
axis at the rear of each plot. The forwardmost line perpendicular to the 
thruster axis represents the plane of the anode pole piece. Measurements 
were made with a moveable Langmuir probe and the recording system 
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described in Reference 6. The raw data were analyzed using the method 
of Strickfaden and Geiler 7 . 

Comparison of the data of Figures 6, 7 and 8 suggests the following 

1. Substitution of higher perveance dished grids for flat grids 

results in: 

a. a reduction in Maxwellian electron density. This re- 
duction was greatest at the 440 ma and least at the 
220 ma flow conditions and most pronounced at each 
flow condition near the screen grid, 

b. an increase in primary electron density. The increase 
was again greatest at the 440 ma flow condition and 
most pronounced at each flow near the screen gird* and 

c. no variation in plasma potential or Maxwellian electron 
temperature that is considered significant. 

2. Throttling the flow into the thruster results in: 

a. a proportionate decrease in Maxwellian and primary elec- 
tron densities throughout the discharge chamber. 

b. a more pronounces plasma potential ridge at the anode, 
and 

c. an increase in electron temperature near the screen grid 
which becomes significant at the lowest flows considered. 

Conclusions 

Installation of high perveance dished grids makes it possible to 
operate a 15 cm diameter SERT 1 1 -type thruster at an ion beam current 
over 600 ma at 85% utilization and 250 ev/ion discharge power with a 
1000 v net accelerating voltage. When the flow rate is throttled to 
obtain a 130 ma ion beam with the dished grids, propellant utilization 
and discharge power change to 60% and 300 ev/ion respectively. This 
degradation in performance with throttling can be reduced if the screen 
grid transparency is reduced with flow rate. Additional testing is 
required to determine if this variation in grid transparency can be 
achieved through variations in grid separation alone. 
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PULSED IGNITION STARTUP STUDIES 


Introduction 

Appreciable cost or weight savings can be achieved in synchronous 
orbit satellite systems if ion thrusters are used to produce the associ- 
ated North-South stationkeeping requirement instead of a chemical pro- 
pellant thruster 8 . For a typical synchronous communication satellite 
having a mass of 700 kg, the 46 m/sec annual North-South stationkeeping 
velocity increment can be met by a mercury ion thruster operating at a 
600 ma beam current level and 1000 volt thruster potential on a duty 
cycle of 230 hours per year. Assuming 60% overall electrical efficiency, 
the thruster could satisfy this requirement by operating for 230 one hour 
intervals from a 1 kilowatt-hour battery which could be charged continu- 
ously from a small solar array (~ 30 watt output). For such a thruster 
system to perform effectively however it must be capable of repeated 
startup which is both rapid and reliable and be able to produce a 600 ma 
ion beam at 1000 volts thruster potential with acceptable propellant 
utilization and electrical efficiencies. Electron bombardment ion 
thruster startup has required resistance heating of the hollow cathode 
tip to a sufficiently high temperature to permit spontaneous ignition of 
the keeper discharge. This mode of startup is frequently unreliable be- 
cause the cathode tip temperatures and the times which are required to 
start the discharge tend to vary with the history of the cathode. Appli- 
cation of a high voltage discharge system developed by Wintucky 9 to 
ignite the keeper discharge is proposed as a means of correcting this 
reliability deficiency. 

Apparatus 

The 15 cm modified SERT II thruster with dished grids, shown in 
Figure 1 and discussed previously, can produce the 600 ma ion beam 
current at 1000 v net accelerating potential necessary to execute the 
North-South stationkeeping mission, and it has been used for this series 
of tests. The high voltage ignition system developed by Wintucky 9 and 
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used in these tests to demonstrate a rapid, reliable startup capability 
is shown in Figure 9. This system utilizes a 1 yf capacitor which dis- 
charges through the secondary winding of a standard 12 v automobile 
ignition coil to produce a high voltage pulse. This pulse is applied 
to a tickler electrode located between the cathode and keeper electrodes 
across which the discharge is to be ignited. Actuation of the pulse 
occurs when the push button circuit is closed thereby triggering the 
silicon controlled rectifier (SCR) C30D shown in Figure 9. 

Oscillographs of the potential difference between the tickler and 
the cathode show the tickler potential rises to voltages in the range 
5000 v tc 10,CCO v in 10 to 30 nsec before breakdown occurs between 
these electrodes. Oscillographs also show multiple breakdown-extinction 
cycles may occur after a single actuation of the circuit and that up to 
several milliseconds may expire between the time of actuation and the 
time when keeper currents in the range of tenths of an ampere begin to 
fl ow. 

The tickler, a 1 mm diameter, 1% thoriated tungsten welding elec- 
trode with its tip ground to a sharp point, was installed perpendicular 
to the cathode axis. The axis of the tickler was located 0.05 cm from 
the face of the cathode tip on the basis of preliminary experiments by 
Wintucky. With the tickler in this location no change in thruster per- 
formance (as shown in Figure 2) was observed. The discharge capacitor 
charging voltage (225 v.d.c.) was selected on the basis of preliminary 
studies on a cathode-keeper assembly which showed this voltage was suffi- 
cient to effect repeatable discharge between the tickler and cathode at 
normal cathode flow rates. The keeper potential of 200 v was determined 
to be sufficient to effect a keeper discharge after the tickler electrode 
had been triggered. A similar preliminary study of the effect of cathode 
tip temperature on keeper discharge ignition probability showed a high 
probability that the keeper discharge would be an arc discharge above 
cathode temperatures around 400°C. Below 400°C keeper discharge igni- 
tions which did occur resulted in a keeper-cathode glow discharge which 
did not result in sufficient cathode heating to bring up the cathode 
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temperature to the point where a transition to an arc discharge could 
occur. 

The only change in the circuity shown in Figure 9 from that sug- 
gested by Wintucky is the addition of the 1000 Q resistor shown in the 
power line supplying the discharge capacitor. This resistor was found 
to be necessary to limit the rate at which the discharge capacitor was 
recharged after the SCR had been tripped. Without the resistor the 
capacitor would frequently not discharge completely for a sufficiently 
long period of time to allow the SCR to reset and in this condition the 
SCR represented a permanent short circuit across the power supply. 

Testing of the system under conditions where breakdown did not occur 
at the tickler electrode resulted in high voltage breakdown through the 
bakelite between the primary and secondary terminals on the ignition coil 
and demonstrated the lowest resistance to high voltage breakdown at this 
point in the system. 

Repetitive Startup Cycle Test Results 

In order to demonstrate the suitability of the thruster and ignition 
system for repeated startup-shutdown cycle operation, two tests were con- 
ducted. The first one involved establishing all mercury flow conditions 
and the neutralizer discharge and then performing 1000 startup and shut- 
down cycles of the main discharge. The second series of tests involved 
establishing all flow rate conditions and then performing 1000 cycles of 
neutralizer followed by main discharge startup with subsequent termina- 
tion of both discharges. 

The main cathode used in these tests was first installed in the 
thruster several months before the tests mentioned above. The cathode 
insert had been dipped in chemical R-500 before cathode installation and 
no additional R-500 was added after this initial coating. The cathode 
had been operated for about 130 hours at the time the first test sequence 
began. The neutralizer cathode is a standard SERT II model; it had been 
operated for over 250 hours and had not been supplied with R-500 beyond 
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that applied at the time of fabrication for the SERT II flight program. 

In conducting the first series of tests the main, cathode and 
neutralizer flow rates were established at 700 ma, 30 ma and ~20 ma 
respectively , and the neutralizer keeper discharge was established at 
a current of 0.2 amps. A main keeper voltage was then established at 
200 v and the arc power supply was set to limit the arc current at about 
3.9 amps once the arc discharge started. The cathode tip temperature 
was next raised to 500°C, a value selected as sufficiently far above the 
necessary temperature for startup (400°C) to assure a high startup re- 
liability. The spark ignition system was actuated next, the keeper and 
arc discharges both established themselves and the keeper current was 
then adjusted to 0.3 amps. The high voltage was turned on next to 
establish a beam current; and the beam current^ arc voltage, arc current, 
keeper voltage and cathode tip temperature were recorded. The high 
voltage, keeper and arc power supplies were subsequently turned off and 
when the cathode tip had cooled to 500°C, the spark ignition system was 
actuated again and the sequence was repeated until 1000 cycles had been 
completed. Two typical startup-shutdown cycles showing the times re- 
quired to complete each portion of the cycle are presented as Figure 10. 

It shows that the typical cycle was completed in less than one minute. 

This test sequence demonstrated the startup probabilities shown below: 

Frequency of startup on first actuation - 99.0% 

Frequency of startup on second actuation - 0.8% 

Frequency of startup on third actuation - 0.2% 

The second startup test sequence required the installation of a 
switch which facilitated connection of the high voltage output lines of 
the spark ignition system to either the main tickler or a neutralizer 
tickler. The neutralizer tickler was constructed and oriented perpen- 
dicular to the neutralizer cathode axis in the same way as the main 
tickler. The test sequence required establishment of the main, cathode 
and neutralizer flows at the same values used for the startup test 
sequence described previously. The neutralizer and main keeper potentials 
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were both raised to 200 volts, the neutralizer cathode was heated to a 
temperature corresponding to a 1.9 amp current through the tip heater 
and the main cathode temperature was raised to 500°C. The arc power 
supply was adjusted to limit the arc current to about 4.1 amps. The 
subsequent events, shown in Figure 11, included initiation of the 
neutralizer discharge using spark ignition, adjustment of neutralizer 
keeper current to 0.2 a, spark initiation of the main keeper and 
arc discharges, adjustment of the main keeper current to a 0.3 a, and 
application of the high voltages to the thruster grids. After data had 
been recorded, the high voltage, arc, keeper and neutralizer keeper 
power supplies were turned off and the main cathode was allowed to cool 
to near 500°C before the startup sequence was repeated. One thousand 
such cycles were executed to complete the entire test. The frequency of 
discharge initiation as a function of number of actuations of the spark 
ignition system is as follows: 



Neutralizer 

Discharge 

Main 

Discharge 

Frequency of startup on 
first actuation 

84.7% 

95.0% 

Frequency of startup on 
second actuation 

5.6% 

3.1% 

Frequency of startup with 
three to five actuations 

4.4% 

1.3% 

Frequency of startup where 
more than five actuations were 
required 

5.3% 

0.6% 


When more than five spark ignition actuations were required to start the 
discharge a component or operator malfunction could generally be identi- 
fied. This malfunctions included: 1) improper positioning of the 

switch used to connect the appropriate tickler to the ignition coil out- 
put, 2) low cathode tip temperature or flow rate, and most frequently 
3) failure of the SCR (Figure 9) to reset following the previous ignition. 
This failure of the SCR was a result of improper design, and it was cor- 
rected by the installation of a 1000 ft resistor in the line from the 
power supply used to charge the discharge capacitor. 
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The variations in neutralizer keeper voltage, keeper voltage, im- 
pingement current, beam current and arc voltage during this test are 
shown in Figure 12 as a function of the number of startup cycles exe- 
cuted. These data show some variations in impingement current, keeper 
voltage, neutralizer keeper voltage and beam current. The variations 
are considered to be due to corresponding variations in flow rates which 
were observed to vary as much as ±15 ma for main flow rates and ±5 ma 
for cathode flow rates during the test. 

Figure 13 shows the condition of the cathodes at the beginning and 
conclusion of each of the tests in addition to a photograph of the main 
cathode tickler at the conclusion of the first 1000 startups. To facili- 
tate comparison all cathodes are oriented the same way in each photo- 
graph. Comparison of Figures 13a and 13b reveals some main cathode wear 
at the edge of the orifice plate (identified on Figure 13b). This wear 
appears at the point where the tickler passed over the orifice plate, 
and it appears to have occurred near the point where the ball of metal 
formed on the tickler (Figure 13d). For the second 1000 startup cycles 
of the main cathode no additional wear is observed, however, and the 
tickler which was replaced for this test appeared to be unchanged over 
the duration of the test. 

The neutralizer cathode shown in Figures 13e and 13f before and 
after the test also shows no significant wear. As can be seen however 
the texture of the orifice plate surface does appear to be different. 

This apparent difference is caused by a coating of R-500 which was placed 
on the cathode face after the test series but before the photograph of 
Figure 13f was taken. The neutralizer tickler showed no significant 
wear after 1000 cycles. 

A comparison of the data on startup reliability of the main cathode 
for the two tests suggests the main cathode was more difficult to start 
during the second thousand tests. After the second thousand cycles both 
cathodes were sent out to be photographed and when they were returned 
neither could be started. Both cathodes were subsequently coated with 
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HIGH VOLTAGE +lkv,-0.5kv NEUT. KEEPER CURR.0.2 a 
TOTAL FLOW -730 ma KEEPER CURRENT 0.3 a 

CATHODE FLOW -30 ma ARC CURRENT 4.1a 
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R-500 and they restarted readily. It was after this that the final neu- 
tralizer photograph was found to be out of focus and it was rephoto- 
graphed to obtain Figure 13f. It is believed that exposure to the 
aromatic chemical vapors in the photographic laboratory where the cathodes 
were stored for about a week while they were waiting to be photographed 
after each test series caused a degradation in startup reliability. 

Rapid Startup Study 

In order to prevent excessive propellant loss during a short thruster 
operating cycle it is necessary to be able to achieve rapid startup and 
shutdown of the system. Shutdown is not considered a particularly severe 
problem because the thruster can continue to operate and produce thrust 
after the vaporizers have been shutdown and until the flow rates drop to 
rather low values. Startup on the other hand requires a substantial flow 
rate to facilitate ignition, and propellant is lost during the time 
period when this flow is being established. Figure 14 shows results of a 
preliminary study indicating the variation in thruster parameters which 
might be expected during a 68 minute period in which the thruster was 
started and operated for an hour. Although the thruster flow rates were 
probably high enough for startup after 4 minutes of this test had past, 
two additional minutes were required to raise the cathode tip temperature 
because of an insufficiently high current setting on the tip heater. It 
is noted that within about three minutes after startup the thruster is 
operating at acceptable conditions, but it is not until 12 minutes after 
startup, when flow rates reach more nearly correct values, that perform- 
ance parameters settle to their steady state values. Flow rates were 
controlled manually during this test, and it is anticipated that an 
automatic control system possibly coupled with higher power heaters 
capable of raising vaporizer and cathode temperatures more rapidly will 
facilitate very substantial gains in this area. Additional testing to 
determine propellant utilization over the complete startup-operation- 
shutdown cycle is required. 
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Conclusion 

A 15 cm mercury ion thruster operating with dished grids is capable 
of producing a 600 ma ion beam at a specific impulse of 2690 seconds and 
a beam ion cost of 250 ev/ion. Rapid, reliable ignition of the neutral- 
izer and main discharges can be achieved using a spark ignition system. 
Use of this ignition system over some thousand cycles results in no 
significant degradation of the thruster components. Installation of the 
tickler electrode results in no increase in beam ion cost or decrease in 
propellant utilization efficiency. Additional attention is needed in the 
area of rapidly acting vaporizers and a control system capable of bring- 
ing flow rates to near the desired values in a time period of the order 
of one minute. 
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SINGLE CUSP MAGNETIC FIELD THRUSTER 
(J. R. Beattie) 

Introduction 

A major factor which determines the lifetime of an electrostatic 
ion thruster accelerator grid is the material erosion which occurs as a 
result of ion-neutral charge exchange reactions. This charge exchange 
reaction rate is dependent on the local ion and neutral number densi- 
ties, and it is desirable to produce a uniform density profile for both 
the ions and neutrals in the vicinity of the accelerator grid to effect 
a uniform erosion rate. As a result then of the ion beam uniformity 
the time integrated thrust for a given grid lifetime is maximized. 

A comparison of the ion density profile existing near the screen 
grid and the beam current density existing near the accelerator grid of 
a Kaufman thruster was conducted by Wilbur 10 and revealed that a good 
correlation exists between these profiles (correlation coefficient of 
0.945) . This high correlation suggests that a uniform ion density pro- 
file near the screen grid should result in a uniform current density in 
the ion beam. 

The only thruster discharge chamber modifications which have re- 
sulted in a significant improvement in the uniformity of the ion beam 
current density profile have been related to the magnetic field configu- 
ration existing in the discharge chamber. The most significant improve- 
ment reported in the literature was achieved by the radial field 
thruster 11 . In this thruster the flatness parameter e (defined as the 
ratio of average to maximum beam current density) was increased from a 
value of about 0.4 for the SERT II thruster to 0.67. The increase in 
the flatness parameter is apparently due to a uniform radial plasma 
density which is found to exist in the radial field thruster. 

Figure 15 presents the critical magnetic field line 12 geometry 
which exists in a SERT II-type divergent field thruster. It is apparent 
from this figure that there should be less ionization occurring at the 
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outer radii of the discharge chamber. This occurs not only because of 
the lower electron densities and energies observed there, but also be- 
cause of the reduced residence time of a neutral mercury atom in the 
ionizing region between the critical field line and the extraction 
grids. Using the residence time approach, one can argue that a properly 
shaped magnetic field would insure that the residence time for a neutral 
atom in the primary electron region (between the critical field line and 
the grids) should be much larger than the time between ionizing colli- 
sions. If t denotes the neutral residence time and v denotes the pri- 
n ^ 

mary electron-neutral atom ionizing collision frequency, the design 
criterion can be expressed as 

t » ~r~ (1 ) 


The residence time is related to the distance d between the critical 
field line and the extraction grid and the average neutral particle 
speed 7 by the relation 


The collision frequency is related to the plasma properties by the 
relation 

v c = n o a o V p » (3) 

where n Q is the neutral particle number density, a* is the electron-atom 
ionization cross section, and V p is the speed of the primary electron. 
Using these expressions for T n and v c the residence time design criterion 
can be related to a minimum distance which should exist between the pri- 
mary field line and the extraction grid. That is 


d 


» X 


o v 


0 ^ 

p 


( 4 ) 


where X + = 1/n o + is the mean free path for ionization of a neutral atom 
ooo 

by a primary electron. Representative values of these parameters near 
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the walls of a SERT II thruster are: n 0 = 10 18 m“ 3 , ? Q = 230 m/sec, 

V p = 3 x TO 6 m/sec (25 ev) , and a* = 2.4 x l(f 20 m 2 . Using these values 

the residence time criterion becomes 

d >> 3.2 mm . 

In order to insure a near 100% probability for neutral ionization 
at the outer radius of the thruster, the value of d should be at least 
an order of magnitude greater than the minimum value given above. That 
is, d should be approximately 32 mm. Since the value of d is estimated 

to be 7.6 mm for the SERT II thruster, the probability for ionization of 

a neutral in this region will be significantly reduced. 

One method of increasing the distance d is by the addition of a 
cusped field electromagnet and anode as suggested in Figure 16. On the 
basis of the increased neutral residence time at the outer radii of the 
discharge chamber, a more uniform plasma density in the vicinity of the 
grids would be expected to exist with the cusped field configuration. 

This would then result in a more uniform ion beam profile as well as a 
reduction in neutral propellant loss at the outer radii of the thruster 
and hence a higher propellant utilization efficiency. An additional 
advantage of this type of field configuration is improved plasma confine- 
ment due to the predominantly axial magnetic field existing in the region 
adjacent to the anode. This improvement in plasma confinement should be 
accompanied by a reduction in discharge power losses plus a reduction in 
radial plasma density gradients provided the magnetic field is confined 
to the plasma region in close proximity to the anode. 

Apparatus 

In an attempt to increase the residence time for neutrals at the 
outer radii of the discharge chamber, the magnetic field modification 
illustrated in Figure 16 was incorporated into a SERT II thruster. The 
major changes to the SERT II configuration consisted of the addition of 
the center magnet pole piece, a reduction in anode length, and the in- 
stallation of the porous screen rear anode. Additional modifications 
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i ncluded the replacement of the SERI II permanent magnets with electro- 
magnets and, during the course of the experiments, the installation of 
a variable magnetic baffle and a high perveance dished grid system. 

The SERT II propellant feed system was replaced by separate main and 
hollow cathode vaporizers; a modification which permitted individual 
control of main and cathode propellant flow rates. A hollow cathode 
with a 0.05 cm diameter orifice was also installed to facilitate sus- 
tained thruster operation at the high ion beam currents anticipated with 
the dished grids. Flat grids were maintained at +3 kv and -1.5 kv and 
dished grids were maintained at +1 kv and -0.5 kv during operation. A 
Faraday cup probe was installed 5 cm downstream of the accelerator grid 
to permit measurement of the ion beam current density profile. In addi- 
tion, moveable Langmuir probes 2 were installed in the main and cathode 
discharge regions to obtain plasma diagnostic data. 

All data presented in this section were obtained at a keeper current 
of 0.3 a. The 37 volt arc voltage operating condition is designated on 
performance curves presented herein by solid symbols. Discharge power 
calculations include keeper discharge losses in addition to the arc 
power losses. 

Early in the test program it was observed that the electromagnets, 
illustrated in Figure 16, were saturating at a relatively low value of 
magnetic field strength. Subsequently the electromagnets were rede- 
signed to increase the cross sectional area off their iron core from 
0.174 cm 2 to 0.658 cm 2 so higher magnetic field strengths could be 
achieved in the discharge chamber. Typical magnetic field measurements 
taken at one point in the discharge chamber with both sets of electro- 
magents operating are presented in Figure 17. These data indicate that 
at a front magnet current of 10 A the magnetic field strength was doubled 
by the use of larger diameter electromagnets and no evidence of satura- 
tion was present. 

An iron filings map of the magnetic field existing in the discharge 
chamber is presented in Figure 18. Magnetic field strength measurements 
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were made at various locations in the discharge chamber for magnet 
currents typical of optimum thruster operation. The results of these 
measurements are illustrated in Figure 19, which shows the magnitude 
(in Gauss) as well as the direction of the magnetic field. The data 
presented in both of these figures were obtained without the magnetic 
baffle energized. 

Anode Configuration Experiments 

One of the first investigations conducted on the cusped magnetic 
field (CMF) thruster was to determine the optimum anode configuration. 
Figure 20 presents performance data obtained at a 312 ma total flow 
rate with flat SERT II grids. Data are presented for the front and rear 
anodes connected electrically in parallel, for the rear anode at float- 
ing potential, and for the front anode at floating potential. These data 
indicate that the best thruster performance is achieved with the parallel 
anode configuration. Also, the qualitative effect of drawing current to 
the rear anode is to reduce significantly the slope of the performance 
curve at high propellant utilization efficiencies. The fraction of the 
arc current drawn to the rear anode also determines the ease with which 
100% propellant utilization is approached. The suspiciously high utili- 
zation indicated in these figures is considered to be a characteristic 
of the CMF thruster design rather than errors in flow rate measurements 
and will be discussed in more detail in a later section of this report. 
Figure 21 presents additional performance data obtained at a high total 
flow rate and with uncompensated high perveance dished grids rather than 
the flat grids used for the data of Figure 20. Data are presented for 
the parallel anode configuration and also for the rear anode at thruster 
body potential. Again it is apparent that the effect of drawing current 
to the rear anode is to reduce the slope of the performance curve at 
high utilization as well as to reduce the baseline discharge losses. 

That the achievement of high propellant utilization is associated with 
the rear anode is clearly demonstrated by the fact that a propellant 
utilization of over 100% was easily obtained with the parallel anode 
configuration . 
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Comparing the parallel anode data at the low total flow (flat grids) 
with that for the high total flow (dished grids) indicates that essen- 
tially the same performance is achieved in either case. This tends to 
support the conclusion reached in Reference 13 and the first section 
of this report that optimum thruster performance over a throttling range 
can be achieved by reducing the grid system transparency as flow rate 
is reduced. 

It was found that the parallel anode configuration resulted in the 
flattest ion beam profile. A comparison of measured ion beam profiles 
for this configuration at high and low total flow rates is presented in 
Figure 22. Comparison of these data with those of Figures 3 and 4 in 
the first section of the report suggests a slightly flatter ion beam 
profile for the CMF thruster. 

Throttling Experiments 

The effects of varying the total propellant flow rate into the dis- 
charge chamber are illustrated in Figure 23. These data were obtained 
with dished grids at a +1.0 Kv screen grid potential and a -0.5 Kv 
accelerator grid potential. The cathode flow rate was in the range 80 
to 110 ma. As the figure indicates, the baseline discharge losses are 
not significantly affected by the total flow rate. However, it is 
evident that the effect of increasing the total propellant flow is to 
flatten out the performance curve and to facilitate operation at propel- 
lant utilization efficiencies in excess of 100%. 

The effect of total propellant flow rate on discharge losses evalu- 
ated at 37 v arc voltage and 85% propellant utilization efficiency is 
illustrated in Figure 24. This curve suggests that the discharge losses 
for these conditions are minimized when the total propellant flow rate 
is in the 550 to 700 ma range. Flow rates in excess of 765 ma were 
found to result in unstable thruster operation. 
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Variable Magnetic Baffle Effects 

A 3.24 cm diameter variable magnetic baffle 14 was installed in the 
CMF thruster to facilitate experimental evaluation of various thruster 
operating conditions. The addition of the variable magnetic baffle per- 
mitted the direct control of the impedance to electron flow across the 
baffle aperture. 

Figure 25 presents a performance comparison between the variable 
magnetic baffle and a 3.53 cm diameter non-magnetic baffle. This larger 
baffle had been used earlier in the Drogram to achieve a 37 v arc volt- 
age at a reasonable cathode flow rate. The figure indicates that the 
performance was essentially unchanged by the addition of the magnetic 
baffle. 

The variable magnetic baffle also permits the acquisition of per- 
formance data over a range of propellant utilization at a constant arc 
voltage. This is accomplished simply by varying the baffle magnet cur- 
rent to maintain a fixed arc voltage as the arc current is varied. 

Figure 26 illustrates the effects on performance of maintaining a con- 
stant 37 v arc voltage throughout the propellant utilization range. 

These data indicate a significant difference in the discharge power 
losses for the variable and fixed arc voltage mode of data acquisition. 

A similar comparison at the higher total flow rates does not reveal much 
difference in discharge losses. This is due to the fact that at the 
higher total flow conditions the arc voltage does not vary significantly 
over a wide range of propellant utilization. 

Cathode Flow Rate Effects 

Figure 27 presents performance data obtained at various cathode 
flow rates and with magnet currents which resulted in minimum discharge 
power losses. These data and those of Figure 25 as well as some data 
obtained without the magnetic baffle and not presented here indicate 
that thruster performance is relatively insensitive to cathode flow 
rate over a wide range of propellant utilization. It has been observed 
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that thruster operation becomes somewhat unstable at very low cathode 
flow rates (-18 ma) and also at very high flow rates {>110 ma) . Based 
on these observations it is concluded that a reasonable cathode flow 
rate range over which thruster performance doesn't vary greatly is 30 
to 100 ma. 

Performance Comparison with SERT II Discharge Chamber 

A performance comparison of the C.S.U. modified SERT II and CMF 
thrusters (both equipped with dished grids) is presented in Figure 28. 
These data indicate that the baseline discharge losses of the CMF 
thruster exceed those of SERT II by about 100 ev/ion at 60% utilization. 
In an attempt to understand the reason for increased losses in the CMF 
thruster a comparison of the measured plasma properties for both 
thrusters was made. The plasma properties obtained from Langmuir probe 
data were averaged along the thruster centerlines and the results are 
presented in Table I. These data indicate that at essentially the same 
flow rate and arc power the ion beam produced by the CMF thruster was 

Table I, Comparison of CMF and SERT II Plasma Properties 


CMF Thruster SERT II Thruster 


Total flow rate 

730 ma 

734 ma 

Cathode flow rate 

85 ma 

34 ma 

Arc current 

3.9 A 

4.32 A 

Arc vol tage 

37 v 

36 v 

Beam current 

449 ma 

632 ma 

Primary electron density 

6.7 x 10 10 cm" 3 

8.1 x 10 10 cm" 3 

Maxwellian electron density 

34 x 10 10 cm" 3 

40 x 10 10 cm" 3 

Maxwellian electron temperature 

2.2 ev 

3.7 ev 

Primary electron energy 

16 ev 

27 ev 

Plasma potential 

38.7 v 

36.3 v 

only about 70% of the SERT II value. 

The reason for 

the significantly 


lower propellant utilization is apparent when one compares the electron 
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densities and energies existing in the two thrusters. From the data 
in the table it can be seen that the CMF thruster primary electron 
energy and Maxwellian electron temperature are about 60% of the SERT II 
values. Also, the electron densities are only about 85% of those 
observed in the SERT II thruster. 

Another striking dissimilarity in the two curves presented in 
Figure 28 is the lack of a "knee" in the performance data presented for 
the CMF thruster. The lack of the "knee" and the apparent propellant 
utilization of over 100% are thought to be intimately related and will 
be discussed in greater detail in the next section. 

Discharge Chamber Theory 

When operating the CMF thruster it was found that the electromagnet 
currents which resulted in minimum discharge losses forced about 70 to 
80% of the arc current to the rear anode. The somewhat unusual perform- 
ance characteristics of the thruster are thought to be due to this large 
rear anode current fraction. 

Consideration of the magnetic field configuration illustrated in 
Figure 18 suggests that a magnetic bottle exists in the region between 
the cathode and center magnet pole pieces. This bottle tends to trap 
the primary electrons which leave the baffle aperture as well as some 
of the Maxwellian electrons which originate from the ionizing collisions 
which occur inside the confinement region. Due to this electron contain- 
ment a high concentration of energetic primaries and Maxwellians would 
be expected to exist in this region. Since the Debye length is very 
small the concentration of ions should be very high also. Physically 
this implies that as the neutral atoms drift through the porous rear 
anode they enter the high energy density region confined by the bottle 
and many become ionized. This argument suggests that the discharge 
chamber might be modeled by a cylindrical region with a planar ion 
source at one end and a set of extraction grids at the other. The mag- 
netic field which exists in the chamber is predominantly axial and tends 
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to inhibit electron flow to the anode. This model of the discharge 
chamber is depicted in Figure 29. 

After the formation of the ion at the source, there are three pos- 
sible methods of extinction. These are: (1) the ion drifts to the 

screen electrode where it either recombines or is extracted, (2) the 
ion recombines on the anode surface, or (3) the ion becomes a double 
ion as the result of an inelastic collision with an electron. Due to 
the rather long distance between the ion source and the screen electrode, 
it appears that the probability of the second and third events occurring 
is quite high. That is, there is a qood chance that the ion formed at 
the source will either be lost to the anode or become doubly ionized. 

Due to the high probability that an ion will be lost to the anode, 
one would expect the discharge losses to be rather high. This is indeed 
the case with the CMF thruster. As was stated above, the effect of the 
magnetic field in the region between the cathode and center magnet pole 
pieces is to confine the electrons. Some electrons will be scattered 
out of this region and it seems reasonable that the electron energy and 
density existing in the remainder of the discharge chamber would be 
relatively low. As a result, the probability of a neutral atom becoming 

ionized in this region of the chamber is reduced. The data presented in 
Table I indicate that the electron energies and densities averaged along 
the CMF thruster centerline are somewhat less than those existing in a 
conventional thruster. Also, for approximately the same input power the 
propellant utilization in the CMF thruster was only about 70 % of that in 
the SERT II. 

Although the region confined by the magnetic bottle has not been 
probed, plasma measurements made elsewhere in the discharge chamber tend 
to support the theory outlined above. Figure 30 presents a comparison 
between the electron number density measured in the CMF and SERT II 
thrusters. The CMF data indicate a general trend for the electron densi- 
ties to increase in the direction of the rear anode. It is also evident 
that the magnitude of the plasma density is somewhat less than that of 
the SERT II throughout the rest of the discharge chamber. 
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Ion Beam Profile Theory 

Since many of the ions formed at the planar source mentioned pre- 
viously are lost to the anode one would expect the ion number density 
to be reduced at the outer periphery of the discharge chamber. This 
effect is confirmed by Langmuir probe measurements taken inside the 
chamber and Faraday probe measurements of the ion beam profile. 

Another factor which can influence the beam profile significantly 
is the presence of double ions in the discharge chamber. Due to the 
higher concentration of ions and electrons in the center of the discharge 
chamber, doubly charged ions are preferentially produced in this region. 
This distorts the single ion current densities at the center of the dis- 
charge chamber and decreases the ion beam flatness parameter. Since 
double ions are produced predominantly by a second ionization of singly 
charged ions this effect is more pronounced in the CMF thruster. In 
this thruster a single ion produced at the rear of the thruster must 
traverse a significant distance before it reaches the extraction grids. 
Hence there is a significant probability that the second ionization can 
occur during this traversal. 

The Effect of Double Ions on Performance Data 

The flat performance curves observed in the cusped field thruster 
require some justification. The following theoretical development is 
intended to show that small double ion fractions result in an apparent 
increase in propellant utilization and an apparent reduction in dis- 
charge power (ev/ion) which causes the performance curve based on 
measured data to be flattened out. 

The ion beam current for the case where both doubly and singly 
ionized species are present is given by the expression 

h = 9 n + v + A o + 2 ^++ v ++ A 0 ’ (5) 

+ t ++ 

where a is the electronic charge, ng and v^ are the number density and 
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velocity of specie £, and A Q in the open area of the grids. In the 
expression above the number densities are taken as the average value 
existing at the accelerator electrode. If there are no double ions 
present the current is given by 

! b+ ■ n/A » {6) 


where n is the average ion number density when there are no double ions 
present. Since double ions are generally created from single ions, the 
specie conservation equation can be written as 


n 


n + + " ++ 


(7) 


Defining the beam current correction factor a as the ratio of the beam 
current when both species are present to the single ion beam current, 
one obtains 


I B +j++ 1 + 2 71 n ++ /n + 

= I B+ ' 1 + n ++ /n + 


(8) 


Where the relation v ++ = J£ v + , which describes the relative velocity of 
double and single ions in the beam, has been used to simplify the results. 
Note that a is given in terms of the ratio of double to single ions 
existing in the beam. In order to relate this ratio to the density ratio 
existing in the discharge chamber one must consider the continuity equa- 
tion. Between two extraction grids one can neglect single and double ion 
production and the conservation of mass (or charge) of a given specie may 
be expressed as 

+ v(nv) = 0 . (9) 

For one-dimensional steady flow this reduces to 

nv = constant . (10) 

Applying the continuity equation to the single and double ions results in 
the relationship 
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n* v* 

M ++ V ++ 

n*v* 


+ ’+ 


n ++ v ++ 

n + v + 


(ID 


where starred symbols refer to discharge chamber conditions at the screen 
grid and the remaining symbols to the beam conditions. The development 
of Appendix A shows that in the discharge chamber v* + - & v*. Using 
this result and the relation v ++ = /l v + the expression above becomes 


n 


★ 

++ 




( 12 ) 


and this can be combined with Equation 8 to yield an equation for a in 
terms of discharge chamber double ion fraction. Figure 31 presents a 
plot of the beam current correction factor as a function of the ratio of 
double to single ions existing in the discharge chamber. 

Since the measured beam current is used in the calculation of the 
discharge losses and the propellant utilization, both of these values 
will be affected by the presence of double ions. The discharge losses 
and propellant utilization calculated from the measured beam current can 
be corrected for the presence of double ions by the relation 

P D = aP D , (13) 

Corrected 

and 

n =-r\ . (14) 

corrected 

From Figure 31 it is apparent that even relatively modest double ion 
fractions can significantly affect the calculated performance. For 
instance, if the beam contained 6 % double ions the effect would be to 
increase the utilization and decrease the discharge losses by about 10%. 
Clearly, the presence of even small double ion fractions has the effect 
of flattening out the performance curve and the effect is more pro- 
nounced at higher propellant utilization efficiencies ( i . e . , larger 
double ion fractions). 
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The Effects of Reducing the Thruster Length 

In Reference 2 it is shown that the baseline discharge losses are 
directly proportional to the ratio of the total area of the primary 
electron region to the open area of the screen grid. For a cylindrical 
primary electron region this ratio is given by 


total 

a 

open 


2ttR 2 + 2uRd 
ttR^Q 


(15) 


where R is the thruster radius, d is the length of the primary electron 
region, and Q is the screen grid open area fraction. For the CMF 
thruster a 50% reduction in the length d would reduce the area ratio 
about 27%. One could then expect a similar reduction in the baseline 
discharge losses. Decreasing d also reduces the magnetic field fringing 
between the front and center magnet pole pieces which might be expected 
to improve the radial uniformity of the plasma density thereby improving 
ion beam flatness. 


It is shown on page 72 that the production of double ions by both 
primary and Maxwellian electrons is a strong function of the volume to 
surface area ratio of the primary electron region. That this is an 
important parameter can be explained on the basis that the production 
of ions is proportional to the volume of the region while the loss of 
ions is related to the surface area confining this volume. For a 
cylindrical region this ratio becomes 

V _ mR 2 d 
A total " 2 " R2 + 


For the CMF thruster a 50% reduction in d would reduce the volume to 
surface area ratio about 30%. A significant reduction in the presence 
of double ions would then be expected as a result. 


Conclusions 

The addition of a center magnet pole piece and a porous screen 
rear anode to a standard SERT II discharge chamber has produced a 
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thruster whose operation is very stable over a wide range of propellant 
utilization. The thruster can be throttled at least over a 2 to 1 range 
in total propellant flow rate at a cost of 50 ev/ion at 85% utilization. 
However, the ion beam current density profile of this thruster has been 
improved only slightly and the baseline discharge losses are signifi- 
cantly higher than those of the SERT II thruster. The typical high flow 
rate performance curve for this thruster design is essentially linear 
and propellant utilization efficiencies of over 100% are easily obtained. 

As a result of the experimental studies conducted on the CMF 
thruster, a fairly comprehensive discharge chamber theory has evolved. 

The theory seems to explain adequately the performance characteristics 
of the thruster in terms of the substantial probability for ions produced 
in the magnetic bottle formed at the rear of the discharge chamber to 
either be lost to the anode or to become doubly ionized. This high 
probability is directly related to the length of the discharge chamber 
and the theory developed would predict that shortening the distance be- 
tween the ion source and extraction grids would reduce both the produc- 
tion of double ions and ion recombination on the anode surface. This in 
turn should result in reduced discharge losses and a more uniform ion 
beam profile. 
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THE PENNING IONIZATION MECHANISM IN AN 
ELECTRON-BOMBARDMENT ION THRUSTER 

(J. Graf) 

Introduction 

Large increases in the ionization rate of a given atomic specie 
can be effected through the addition of a second specie which has a 
metastable level above the ionization potential of the first. This 
phenomenon is known as the Penning effect, and ionization is effected 
by a metastable specie two atom which collides with a specie one atom 
to produce a specie one ion and a ground state specie two atom. Since 
argon and helium have metastable levels above the ionization potential 
of mercury a collision between either argon or helium metastables and 
a mercury atom should theoretically result in a mercury ion. Since 
the cross section for this reaction is of order 10 14 cm 2 , which is 
two orders of magnitude above the cross section for mercury ion produc- 
tion by electron bombardment, it was hoped the Penning interaction would 
supplement mercury ion production by electron bombardment and hence in- 
crease the ion production rate significantly. 

Theory 

During the Penning interaction in a mercury bombardment thruster, 
an excited helium or argon atom for example transfers its energy to the 
mercury atom. The amount of energy transferred (19.8 ev for helium and 
11.53 ev for argon) is higher than the first ionization potential of 
the mercury (10.43 ev) and so ionization is possible. For the case of 
helium which will be used as an example here because of the availability 
of cross section data the process is described by 

He 111 + Hg ■+ He + Hg + e . (17) 

The metastable state is one which does not readily radiate to 
ground level via emission of a quanta of radiation and it therefore has 
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a long radiative lifetime (on the order of 10 -1 sec). If the time for 
Penning interaction were much less, then this would be the dominant 
de-exitation mode. For helium this interaction time is defined by the 
equation 


1 v He n Hg ov He 

where A p enning is the mean free path for the Penning interaction, v He 
is the helium velocity, is the mercury density and a is the Penning 
cross section. Using the average thermal velocity of helium at 50CTK, 
a Penning cross section of 1.4 x 10 14 cm 2 and a mercury density of 
!0 18 atom/m 3 the interaction time, t. is 4.4 x 10 4 sec. Thus, during 
its existence, the metastable has a high probability to collide with 
and ionize a mercury neutral. 

To estimate the effect of the metastable helium on the mercury ion 
population, one must first calculate the metastable density. The steady 
state, metastable mass balance for an elemental volume in the plasma is 


Creation 

excit m 

a i n He n e + e ! n He = “ 3 "Hg n He 


Destruction 
+ m 

+ a, * n Hg n He 


+ a s n He "He + a s n He n He 


+ a7 "e "He + n He 


aj = rate coefficient for neutral helium-electron collision 

a 3 = rate coefficient for metastable helium-mercury atom collision 

a 4 = rate coefficient for metastable heli urn-mercury ion collision 

a 5 = rate coefficient for metastable helium-helium atom collision 

a 6 = rate coefficient for metastable helium-helium metastable collision 

a 7 = rate coefficient for metastable helium-electron collision 

3i = rate coefficient for metastable helium production from excited states 

B 2 = rate coefficient for metastable helium-wall collision 
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Most work done on the production of metastable helium has combined 
the cross section of excitation due to direct electron bombardment with 
that due to step-wise metastable production from a higher energy state. 
Thus, rate coefficients a! and can be combined into a third coeffi- 
cient, a 2 s the rate coefficient for total metastable production 

«m He n e + 6i n** clt - ct 2 n He n e • (20) 

(The coefficient 3 X is a function of electron density since the neutral 
helium is excited initially by electron impact.) 

Assuming an infinite plasma, one can neglect de-excitation due to 
the presence of the wall (e 2 = °)- This is not strictly valid as the 
Penning interaction length is on the order of the thruster dimension, 
approximately 0.71 meters/collision. For this preliminary study, how- 
ever, the above assumption will be retained. In order to estimate the 
mercury production rate from metastable interactions an electron and ion 
density of 10 17 m' 3 and a neutral mercury density of 10 18 m“ 3 will be 
assumed. 

Figure 32, taken from Massey and Burhop 15 , shows that the maximum 
cross section for production of 2 3 S helium metastables (19.8 ev) occurs 
at an electron energy of about 40 ev--approximately the energy of a 
primary electron in the thruster chamber. Helium has another metastable 
level, 2 1 S, at 20.6 ev, which, according to Biondi 16 , decays rapidly, 
via electron collision, to the 2 3 S state 

He m (2 1 S) + e' -*• He m (23S) + e" + 0.79 ev . (21) 

Figure 32 also shows that the production of the 2*S metastable is 
not at a maximum at the thruster operating voltage. For these two 
reasons, the assumption will be made that only the 2 3 S (19.8 ev) meta- 
stables need be considered. The maximum numerical experimental cross 
section (a 2 ) associated with Figure 32 was determined initially by 
Maier-Leibnitz 17 and found to be 5 x 10“ 18 cm 2 . Other experimenters 
found the cross section to be 4(±.3) x 10“ 18 cm 2 (Schultz and Fox 18 ), 
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2.6(+0.4) x 10~ 18 cm 2 (Fleming and Higginson 19 ) and 3(±0.7) x 10 18 cm 2 
(Holt and Krotkov 20 ). The order of magnitude agreement is sufficient 
for this rough analysis. 


The neutral helium density in the discharge chamber, n Hg can be 
calculated from simple kinetic theory, assuming helium is lost through 
the thruster grid system at the supply rate and that the loss rate is 
determined by the sharp edged orifice expression for free molecular 


fl ow 


4m. 


n 


He 


He 7 A 
o 


( 22 ) 


m u , the helium atom flow rate is assumed to be 100 ma equivalent 

n6 * 

(6.25 x 10 17 atom/sec) and A, the open area of a SERT II grid system, 
is .0124 m 2 . the average velocity of the neutral helium, v Q , is given 
by the kinetic theory expression which at the thruster temperature of 
500°K has the indicated numerical value. 

v 1 .63 x 10 3 m/sec (23) 

o yj irm He 

Substituting these values into Equation 22, the neutral helium density 
is determined to be 1.24 x 10 17 m 3 . 

The rate coefficients (a) are defined to be the product of the 
associated cross section and the relative velocity of two colliding 
species, v r , which for the case of a colliding electron and heavy par- 
ticle is essentially the electron velocity. The relative velocity be- 
tween two heavy particles can be determined from the equation 21 



where m* is the reduced mass. It has been shown previously that the 
maximum production of 2 3 $ metastables occurs at about 40 ev. Thus, for 
this rough analysis the electron energy will be assumed to be 40 ev 
where the electron velocity is 3.75 x 10 6 m/sec. This corresponds to a 
rate coefficient for metastable production, a 2 , of 1.875 x 10 15 m 3 /sec. 
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The right hand or destruction side of Equation 19 is much more com- 
plicated. Metastable helium can be destroyed by collisions with mercury 
neutrals and ions, helium neutrals and metastables and electrons in 
accordance with the following equations: 


He m + Hg -*• He + Hg + + e 

(25) 

He m + Hg + -*• He + Hg ++ + e 

(26) 

He m + He 2He + e 

(27) 

He m + He m ■» He + He + + e" 

(28) 

He m + e" -* He + e" "l 


He m + e' He + + 2e~J 

(29) 


The cross section for the Penning effect coefficient, a 3 , has been 
determined by Biondi 16 to be 1.4(±0.3) x 10“ 11+ cm 2 . The author has 
found no mention of process 26 in the literature; theoretically how- 
ever the ,19.8 ev metastable should be capable of ionizing the mercury 
a second time (18.75 ev ionization potential). Thus, for this rough 
analysis, the cross section will be assumed to be that or process 25, 

(a 4 - 1.4(±0.3) x ICf 1 ' cm 2 ). 

The cross section for helium neutral and metastable collisions is 
small, partly due to excitation transfer, i.e., a metastable atom trans- 
ferring energy directly to another atom without radiative losses. 

Biondi 16 calculated a value of 9.6(±0.1) x 10~ 21 cm 2 but it is not be- 
lieved that this is the 2 } S metastable neutral cross section. Allison 
et al 22 repeated the calculation for 2 3 S mstastables and found 
a 5 * 2 x 10' 20 cm 2 . 

Much work has been done on the metastable - metastable cross section 
and they agree on the order of magnitude, 10 cm 2 . Specifically, the 
value found by Hurt was a 6 = 1.2 x 10 cm 2 . This value, although of 
the order of a 3 , will not have a large effect in the chamber due to the 
low helium metastable density compared to that of mercury. However, much 
energy could be lost in the region of metastable production due to this 
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process. Geometries which spread out the metastable production region 
may prove to lower this loss if, in fact, it is significant. 

The final process, metastable de-excitation by electron collision 
can, according to Loeb 24 , result in a helium neutral or ion. Bates 
et al 25 calculated this interaction coefficient by inverting the 
Schul tz-Fox 18 curve for the 2 3 $ metastable production. The coefficient 
was found to be a 7 = 8.4 x 10" 16 m 3 /sec. 

Using these cross sections and heavy particle velocities corres- 
ponding to a wall temperature of 500°K one calculates the following 
rate coefficients: 

a 3 = 2.3 x 10” 15 m 3 /sec 

c <4 = 2.3 x 10’ 15 m 3 /sec 

a 5 - 4.6 x 10“ 21 m 3 /sec 

a 6 = 2.75 x 10‘ 15 m 3 /sec 

a 7 - 8.4 x 10‘ 16 m 3 /sec 

Collecting the interaction coefficients and the various densities, 
substituting into Equation 19 and solving the resulting quadratic, the 
helium metastable density was found to be n™ g = 8.81 x 10 15 m 3 . Thus, 
the metastable density is two orders of magnitude below the neutral 
helium density. When combined with the neutral mercury density and the 
Penning coefficient, a 3 , the mercury ion production rate due to helium 
metastables is ct 3 n jj e n Hg = 2.03 x 10 19 m’ 3 sec“ 1 . 

When this rate is compared to the neutral mercury-electron inter- 
action ion rate at 40 ev, it is found to be three orders of magnitude 
lower, namely, 

a 8 n Hg n e = 2.25 x 10 22 m“ 3 sec -1 . 

Remembering that this analysis did not include the interaction of the 
metastables and the wall, one concludes the metastables should have 
little effect on the rate of ion production. For the metastable 
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produced ions to have an effect the Penning production rate must be of 
the order of the mercury-electron production rate. This implies a 
helium mass flow rate that is four orders of magnitude higher (m He = 

1.02 x 10 6 ma). 

Experiment 

To confirm the preceding theory, a Kaufman thruster was modified 
to permit the introduction of the seeding gas through the hollow cathode. 
When helium was used, high cathode flow rates (> 1 amp eq.) were re- 
quired to sustain the keeper discharge. In addition, the discharge was 
found to be unstable and shorting was observed between the grids. (The 
shorting was a result of high test chamber pressure due to the high 
helium flows.) As there was great difficulty in sustaining the arc dis- 
charge, no reliable data could be recorded. 

When argon was used rather than helium, the thruster operation was 
found to be stable. The argon flow rate was also very high (> 1 amp eq.) 
and any attempt to reduce it extinguished the arc. This resulted of 
course in poor utilization efficiency, and the observed baseline dis- 
charge power of 316 ev/ion was greater than that observed with no argon 
f 1 ow . 

Since the above configuration did not provide any satisfactory re- 
sults, a third feed tube was installed to facilitate helium or argon 
distribution with the main mercury flow, throughout the discharge cham- 
ber. The cathode operated on its own independent mercury reservoir and 
this design permitted accurate measurement of all flows. With this 
configuration, stable thruster operation could be established and the 
Penning gas could then be introduced in amounts from trace flow rates 
up to a few amps equivalent. The test results and performance data 
obtained with argon are shown in Table II while those obtained with 
helium are presented in Table III. The utilization efficiency was cal- 
culated using the sum of the Penning gas flow and both mercury flows. 
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Conclusion 

The results displayed in Tables II and III show no improvement in 
the performance characteristics of the thruster when either helium or 
argon are introduced into the discharge chamber. This confirms the 
theoretical prediction that the metastables have no substantial effect 
upon the mercury ion population. 
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Table II 

Mercury-Argon Operating Conditions and Performance Data 


Mercury Flow Rates 
Cathode - 30 ma 
Main -700 ma 



Keeper Current 

0,3 a 

Run Number 

1 

2 

3 

4 

5 

Argon Flow Rate (ma) 

0 

57.8 

422 

>870 

>870 

Arc Current (a) 

4.27 

4.25 

4.26 

4.20 

4.25 

Arc Voltage (v) 

37.7 

38.5 

37.9 

38.3 

38.1 

Keeper Voltage (v) 

11.2 

11.4 

11.3 

11.4 

11.1 

Beam Current (ma) 

609 

609 

614 

610 

628 

Impingement Current (ma) 

2.4 

2.4 

2.4 

2.4 

2.5 

Beam Ion Cost (ev/ion) 

269.8 

274.2 

268.5 

269.3 

263.1 

Propellant Utilization 
(*) 

83 

77 

53 

-- 

— 


Table III 

Mercury-Helium Operation Conditions and Performance Data 

Mercury Flow Rates Keeper Current - 0.3 a 

Cathode - 30 ma 
Main -700 ma 


Run Number 1 2 3 4 5 


Helium Flow Rate (ma) 

0 

82.2 

233 

500 

>10 3 

Arc Current (a) 

4.22 

4.24 

4.24 

4.23 

4.24 

Arc Voltage (v) 

38.8 

38.9 

38.9 

38.9 

38.9 

Keeper Voltage (v) 

11.8 

11.9 

11.9 

11.9 

11.8 

Beam Current (ma) 

609 

616 

615 

615 

617 

Impingement Current (ma) 

2.6 

2.6 

2.6 

2.6 

2.6 

Beam Ion Cost (ev/ion) 

274.7 

273.5 

274.0 

273.4 

273.1 

Propellant Utilization {%) 

83 

76 

64 

50 

— 




DOUBLE ION PRODUCTION IN MERCURY 
ELECTRON-BOMBARDMENT THRUSTERS 


Introduction 

Up until the past few years concentrations of doubly ionized mercury 
in electron bombardment ion thrusters were considered small when primary 
electron energies were below about 35 ev because of the relatively low 

4 2 

cross sections for double ionization of neutral mercury at these energies 
Subsequently investigation of plasma within thrusters have included 
spectroscopic and mass spectrometer measurements which suggest the double 
ion densities are considerably higher at these low energies than analysis 
based on this reaction would suggest 27 ’ 28 . An alternative reaction which 
might account for the higher than anticipated double ion densities is 

Hg + + e” -► Hg ++ + 2e . (30) 

This reaction can be significant in mercury bombardment thrusters if 
the singly charged ion densities are sufficiently high. Bombarding elec- 
trons which could be important in this reaction are the monoenergetic 
primary electrons and those with a Maxwellian energy distribution. If 
one assumes ionization from various excited states is negligible double 
ion concentrations resulting from the above reaction as well as direct 
double ionization of neutral mercury can be calculated and the results 
can be compared with experimental measurements. 

Theory 

Neglecting ionization of excited neutrals and volume recombination 
of double ions s the generation rate of single ions R* in a volume V of 
uniform atomic density n Q , Maxwellian electron density n m and distribu- 
tion function f(E), and primary electron density n p and energy ? is 
given by 

00 

R g ’ % V ["p % <5 p > v p U p ) + f 4,n m a + 0 (E)f(E) c3(E)dc] , (31) 

0 
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where v p =^y is the speed of primary electrons, c(E) = y ~ is the 

speed of Maxwellian electrons at energy E, and cr*(E) is the cross sec- 
tion for single ionization from the neutral ground state by an electron 
of this energy. The first and second terms of Equation 31 corresponds 
respectively to single ion production via primary electron bombardment 
and Maxwellian electron bombardment. The Maxwellian distribution 
function 29 f(E) is given by 

f < E > =(4 t) 3/2 ex p(r) ■ (32) 

where m and T are the electron mass and temperature and k is the Boltzmann 
constant. 

Single ions are assumed lost from the volume only as a result of 
migration through the surface of area A which bounds the volume and 
through ionization to the doubly ionized state. Since double ionization 
is effected again by both primary and Maxwellian electrons the singly 
ionized ions loss rate is given by the equation 

00 

Rj = n + v + A + n + v|n p CT| + ( 5 p )v p ( 5 p ) + j 4Kn m f(E)a+ + (E)c3(E)dcj 
** 0 

( 33 ) 

where n + is the single ion density, a* + (E) is the cross section for 
double ionization from the singly ionized state by an electron of energy E, 
and v + is the single ion velocity across the surface confining the plasma, 
determined by the modified Bohm criterion 30 . The fact that v + is not 
altered because of the presence of the double ions in the plasma is dis- 
cussed in Appendix A. 



In the steady state the production and loss rates are equal and one 
can therefore combine Equations 31 and 33 and solve for the single ion 
to neutral density ratio in terms of measured and calculable thruster 
paramete rs 
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n+ -< 

r n 0 

</(£ ) v (E ) + / 
n m 0 p p p J 

° r + 1 
47Ta Q (E)f (E)c 3 (E)dc 

> 


v. A n 

-t- _ + _EL a (r ) v (r ) + 

Vn m n Q + u p ; W 

^ m m { 

00 

j" 4ira^ + (E)f(E)c 3 (E)(ic^ 
D 


( 35 ) 


In a similar manner one can determine the double ion generation rate R q 
due to both direct ionization of neutrals and ionization of single ions, 
by both primary and Maxwellian electrons. 



"c v {vo + <V VV + / 4 ™A +(E)f(E)c3(E)dc } 

^ 0 

00 

+ n + V | r p a + + (^p) v p<Sp) + / ^ 7ra + + (E)f(E)c 3 (E)dc|. 


(36) 


The loss rate for double ions is assumed to be due only to migration 
across the bounding surface of area A and is given by 


R 


++ 

£ 


n ++ v ++ A 


(37) 


The double ion velocity v ++ necessary to assure a stable sheath is larger 
by the factor JT than that given by Eauation 34 as discussed in Appendix A. 

When Equations 32, 35, 36 and 37 are combined, an expression for the 
ratio of doubly ionized to neutral atoms can be obtained which is a func- 
tion of atom and electron densities, temperatures and energies; the 
volume to surface area ratio of the plasma interaction region; and the 
cross sections for the various reactions. The following equation shows 
how these parameters appear in the final result; for simplicity Equation 
32, the distribution function expression, has not been substituted. 



-72- 


n, , Vn 
++ m 


n Av,, 

o ++ 


oo 

a+ + v + f 4wa* + f c 3 dc 

n m 0 P J ° 


(38) 


n c 

-£■ a + v + / 4ircr f c 3 dc 

"m 0 P J 0 






w~ + % a + + v p + / 4mJ + +f c3dc 


> 


m m 


oo ” 

| ^ c,; + v p + f 4 TO ^ + f C 3 dcj 


J 


This equation can be solved readily on the computer once the cross section 
variation with electron energy is specified. 


The experimental electron bombardment cross sections for the forma- 
tion of single and doubly charged ions from neutral ground state mercury 
are available from References 26, 31 and 32, and those for ionization from 
the singly to doubly ionized state are given in References 33 and 34. 

These cross sections were approximated using a Gryzinski type expression 
and the extent of the agreement between these approximate cross sections 
and the measured ones is presented as Figure 33. 


Resul ts 

Figure 34 presents typical double ion concentrations determined from 
Equation 38 as a function of primary electron energy. The solid line 
represents what is considered a base operating condition with the volume 
to area ratio typical of a 30 cm dia. thruster. The other curves show 
the effect of varying the indicated parameter. Comparison of the dashed 
and solid lines shows for example a somewhat uniform reduction in double 
ion concentration occurs over the primary electron energy range when the 
Maxwellian electron temperature is reduced from 7 ev to 4 ev. The curve 
identified by a centerline corresponds to a 15 cm dia. thruster with 
otherwise similar conditions to those of the solid line. The smaller 
thruster is observed to exhibit much smaller double ion fractions than 
the 30 cm dia. one. This suggests why the presence of double ions in 
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thrusters became much more evident recently when active development of 
larger thrusters began. One observes a variation in primary-to-Max- 
wellian electron density ratio results in a profound effect at high pri- 
mary electron energies and as one might expect a smaller effect at 
energies at and below 20 ev. A 50% increase in Maxwellian and primary 
electron densities is observed to bring about a 100% increase in double 
ion fraction over the entire range of primary electron energies. A re- 
view of these results also shows that the dominant mechanism for produc- 
tion of double ions is the two step process involving ionization of 
singly ionized mercury into the doubly ionized state. 

The concentration of single ions can be calculated from Equation 35. 
Results of this calculation are presented in Figure 35 for the same input 
parameters as Figure 34. The calculated single ion fractions are observed 
to fall into the range generally considered to be reasonable for thruster 
plasmas . 

Comparison with Experimental Results 

Milder and Sovey 27 have measured double ion fractions (n ++ /n o ) using 
a spectroscope on an operating 30 cm diameter ion thruster. Assuming the 
plasma properties suggested in their paper (Maxwellian electron density 
of 10 17 m~ 3 and temperature of 7 ev and primary electron density of 
5 x 10 15 rrf 3 ) and volume to surface area ratio appropriate to the 30 cm 
dia. thruster (0.033 m) their data lend themselves to direct comparison 
with results of this analysis. Figure 36 presents this comparison as a 
function of primary electron energy; the Milder-Sovey thruster center- 
line double ion fractions are shown as solid data symbols confined by 
error bars and the theoretical results are shown as the solid line. The 
significant discrepancy which is evident from this figure would be elim- 
inated if the Maxwellian electron density were 1.5 x 10 17 rrf 3 and the 
primary to Maxwellian electron density ratio were increased to 0.15 as 
indicated by the dotted line. Such an adjustment of variables can be 
argued as necessary to bring the single ion to atomic density ratios 
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( Figure 35) nearer to the 10% values considered typical of ion thruster 
discharge chambers. Additional improvement in the agreement between 
these theoretical results and the data of Reference 27 may be possible 
if ionization routes through metastable and resonance excited atomic 
states are considered. 
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THERMAL FLOW METER 


Introduction 

In the course of conducting experimental studies of mercury bombard- 
ment thrusters, considerable time is consumed in waiting to obtain 
accurate mercury flow rate information. In addition, mercury flow measure- 
ment techniques presently in use involve a measurement of mass or volume 
of mercury, and they therefore provide integrated flow rate information 
which implies an average flow rate. It is very difficult to determine 
flow rate trends from such information. These problems could be alle- 
viated if a continuously indicating, accurate flow meter were available. 

The thermal flow meter described herein, which is based on the concept 
suggested by Laub 35 , does indicate mercury flow rates continuously over 
the range of interest for ion thruster testing with acceptable accuracy. 

Theory of Operation 

The basic concept on which this device operates is illustrated in 
Figure 37. As liquid mercury passes through the flow tube toward the 
vaporizer and thruster, it is heated slightly by the fixed power heater. 
Mercury entering and leaving the flow meter is maintained at ambient 
temperature by a heat sink. The proximity of the high temperature sensor 
to the heater and the fact that the temperatures at opposite ends of the 
tube are held at ambient values implies a temperature difference between 
the two sensors shown in the figure. These sensors (thermistors) are 
connected into a bridge circuit which exhibits an output proportional to 
the tenperature difference which they sense, and this output can be re- 
lated to flow rate. 

A simple one-dimensional theoretical model of this flow meter con- 
cept has been developed 2 , and it has been used to design such a device. 

The flow meter, described previously in Reference 2, consists of a 0.3 mm 
i.d. teflon flow tube fitted with a 60 milliwatt heater and 1.1 mm dia. 
glass bead, 10 Kft thermistors all of which is enclosed in a vacuum iso- 
lation chamber. The flow tube is fitted within a vacuum feed-through 
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in the manner suggested by Figure 38. 

The output from the device is obtained through the basic bridge 
circuit shown as Figure 39. Initial tests were made with 6.8 ICQ 
resistors as the components RJ, ^ in this figure, and it was observed 
that the output of the device drifted with ambient temperature fluctua- 
tions. The reason for this may be understood by considering the follow- 
ing: 

1) The resistance of a thermistor may be assumed to vary with 
temperature in accordance with the equation 

R = R 0 e" aT (39) 


where R Q and a are constants. Hence the resistance R ref at an ambient 
temperature T g is 


R ref ' R o e ' aTa 


(40) 


and the resistance R m= „ some A degrees hotter than T is 

MIS aS G 


R 


meas 


R e” 0l(T a + &) . 


(41) 


2) The output of the bridge circuit is given to first order by 
the equation: 


E 

out 

^battery 


^ref ^ref 

*Yef + R ref ^meas + *Vef 


(42) 


3) Substitution of Equations 40 and 41 into 42 yields an expres- 
sion for output voltage (E QUt ) which is a function of T a for R' rgf 
constant. This implies the observed drift. 

If one substitutes thermistors which sense the ambient temperature 
for the resistors R' f , and if all four thermistors are assumed to be 
identical (R and a identical for each) then the combination of Equations 
40, 41 and 42 yields 
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E out ] 1 

E battery 2 1 + e ' 


(43) 


This equation is independent of T a and hence no drift should be 
observed. Substitution of thermistors into the bridge circuit for the 
components in the experimental system did result in a reduction in 
drift, but it did not eliminate it completely. Presently four matched 
thermistors are being purchased to eliminate possible drift because of 
variation in the constants R 0 and a among the thermistors. 


Experimental Performance 

The flow meter using unmatched thermistors has been calibrated by 
installing it in the line between a 0.9 mm di a . bore glass flow tube 
and the main vaporizer of a thruster and measuring the time required for 
the mercury meniscus to fall between two scribe marks confining a known 
volume of mercury while the flow meter output remained essentially con- 
stant. The resulting calibration curve for the thermal flow meter is 
presented as Figure 40. It shows a nearly constant sensitivity of about 
10 mv/gm/hr over a 0.2 to 5.0 gm/hr flow rate range, and it appears higher 
flow rates could be measured. It is apparent from the curve that the 
flow meter performs well at somewhat higher flow rates than had been con- 
sidered in the initial design. This suggests the conduction through the 
mercury may be higher than had been anticipated and that the somewhat 
higher mercury flow velocities are necessary to achieve significant con- 
vection alterations of the otherwise linear conduction temperature pro- 
file. By increasing the sensitivity of the readout equipment the flow 
meter appears to perform satisfactorily at low flow rates (< 0.3 gm/hr.). 
The scatter in the data points of Figure 40 is considered to be due to 
variations in ambient temperature and scatter in flow rate measurements 
obtained from the flow tube. 

The linear relationship between flow rate and bridge output voltage 
shown in Figure 40 is noteworthy. It occurs because the mercury tempera- 
ture near the flow meter heater is inversely proportional to flow rate 
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while the bridge output (Equation 43) can be set up so, to first order, 
its output is inversely proportional to temperature. These facts com- 
bine to effect a bridge output which is directly proportional to flow 
rate. 

Several months of operating experience with the flow meter have 
demonstrated that it is very useful as a tool for indicating flow rate 
trends. This usefulness is illustrated in Figure 41 which displays flow 
meter output as a function of time for a flow rate initiation - stabiliza- 
tion - termination cycle. After the flow meter heater had been allowed 
to stabilize and the zeroing potentiometer (Figure 39) had been adjusted 
to indicate zero output at the existing no-flow condition, the vaporizer 
was energized. The figure shows a large negative flow rate associated 
with expansion of the mercury occurs immediately after the vaporizer 
heater is energized. This is followed by the establishment of a posi- 
tive flow rate about six minutes after vaporizer energization. After 
a period of vaporizer heater current adjustment which was necessary to 
achieve the desired flow rate, stabilization at a steady flow rate is 
observed (about 15 minutes after vaporizer heatup began). Because of 
the one minute time constant associated with the flow meter, the indi- 
cated flow rate is in good quantitative agreement with the true flow 
rate only in this steady state region. When the vaporizer heater is de- 
energized an increase in flow rate, associated with contraction of the 
mercury as it cools, is observed. Subsequently the output drops to zero 
exponentially as the flow meter temperatures stabilize back to their 
no-flow values. 

Figure 42 displays flow meter output as a function of time while 
the thruster is operating. Two interesting details are apparent from 
this figure. They are: 1) the tiny blip identified with high voltage 

arcing which is observed occasionally and 2) the increase in flow rate 
observed when the level of mercury is raised in the flow tube thereby 
increasing the pressure in the feed system. These are interesting 
because they indicate 1) normal high voltage arcing does not disrupt 
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flow meter operation significantly (no isolator was installed in the 
feed line for these tests) and 2) changes in the mercury supply pres 
sure can alter the associated flow rate. This later phenomenon is 
considered to be a result of vapor at the liquid mercury-vaporizer 
interface which affects the rate of heat transfer to the mercury and 
hence its rate of vaporization and transfer through the porous vapori 
zer plate. 
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APPENDIX A 

ION VELOCITIES TOWARD BOUNDARIES FOR PLASMAS 
CONTAINING DOUBLY CHARGED IONS 


The basic Bohm criterion 36 which determines, in a Maxwellian plasma, 
the ion velocity necessary to assure a stable plasma sheath at a boundary 
has been modified by Masek 30 to account for the presence of a mono- 
energetic electron group in addition to the Maxwellian group. In view 
of the small but significant double ion concentrations being measured 
in ion thrusters it is appropriate to determine if any additional modifi- 
cation is needed to account for this. 

The one-dimensional Poisson equation that applies in the sheath 
region is 

= 4ire(2n ++ (x) + n + (x) - njx) - n p (x)) (Al) 


where V, the negative of the potential in the sheath is taken as zero at 
the boundary, e is the electronic charge, x is the position variable, and 

n . n l5 n and n are the double ion, single ion, Maxwellian electron, 

++ + m p 

and monoenergeti c electron densities respectively. The Maxwellian and 
primary electron densities are given b.y 30 

n m (x) = n m (0) ex p(^) (A2) 

n (x) = n (0) (1 - J-) (A3) 

P K 


where n (0) and n (0) are the electron densities at the plasma-sheath 
m p 

interface, T is the Maxwellian group temperature, k is the Boltzmann 
constant and E p is the monoenergeti c group energy. The singly charged 
ion density, as in the original analysis 36 , is given by the expression 
for conservation of current density 


r E + l 1/2 

n + = n +(°) Le^TT 


(A4) 
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where n + (0) is the ion density at the plasma-sheath interface and E + is 
the ion energy (in ev) at the interface. Similarly the double ion den- 
sity is given by 


n ++ = n ++(°) 


i 1 / 2 


L E ++ 


'++ 

+ V 


(A5) 


where E ++ represents the voltage through which the double ion fell to 
achieve the energy it has at the interface (multiplying E ++ by twice 
the electronic charge yields the double ion energy). 

Substituting Equations A2 through A5 into A1 one obtains 


d 2 V 

dx 2 




- n p(° ) ( 1 -^)], (A6 

Multiplying Equation A6 by dV/dx, integrating and solving for the con- 

dV 

stant of integration to satisfy the boundary condition ^ = 0 at V = 0 
yields 


(f ) 2 - 

+ n m (0) (T)( exp (' If) • 

Near the plasma-sheath interface 



V « « £ 




(A8) 


The terms in Equation A7 pertaining to the ions and Maxwellian electrons 
may be expanded in a series and higher order terms, determined by Equations 
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A8 may be neglected and this equation becomes 



Recalling that charge neutrality exists at the olasma -sheath boundary 
the equation 


2n ++ (0) + n + (0) - nJO) + n p (0) , (A10) 


may be incorporated into Equation A9 to yield 


dV\ 2 . 
,dx/ 

dV 


8tt6 


n ++ (0) n + (0) _ %(°) e n p (°) 

" "2E^ 4E7" + 2k + 2E 


V 2 


(All) 


In order for to be real, one must require 
dx 



n (0) 

n7(oy 


1 \ n ++ (°> 




(A12) 


The second term of Eouation A12 is of order [0.1 e/kT] since E p » e/kT 
and this term may therefore be neglected. For sufficiently small double 
ion fractions the third term is also small and one obtains from Equation 
A12 the required energy of the single ions 


p /kT\/ n +(°) 

E ' M2i7 rTTO] 



n ( 0 ) + n„(0)' 
p' ' m v ' 

n m (0T 


(A13) 


This is the same as the result obtained by Masek for no double ions, and 
it corresponds to an ion velocity of 


v 


+ 



n p (0) + n m (0)^ 

"7 


(A14) 


How small the ratio n ++ (0)/n + (0) must be is determined by the magnitude 
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of E ++ /E + . In order to determine the magnitude of this ratio consider 
n + (0) small compared to n ++ (0) in Equation AH and the criterion 
becomes 


n ++ (0) 


E ++ - n m (0) 
nr 



(A15) 


n (0) + n (0) 

Since n ++ (0) = ^ for charge neutrality at the interface the 

condition becomes 


E 


++ 



n p(°) + "m (0) \ 

"jor — ) 


(A16) 


Equations A13 and A16 taken together suggest E ++ /E + is of order unity and 
therefore that the third term of Equation A12 could be neglected until 
the n ++ (0)/n + (0) reaches a value of the order of 10%. 

Equation A16 which sets a limit on the doubly charged ion velocity 
to assure a stable sheath for the case of negligible single ion density 
yields a value of 


'++ 



= JT v 


+ * 


< 


(A17) 
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